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I. INTRODUCTION 

Electron transport through few electron islands is 
strongly influenced by electron correlations between its 
atomic like states. A variety of correlation effects, in- 
cluding Kondo effec t 1 ' 2 ' 3 and resonant current peaks^ 
have been observed. For low conductance of the tunnel 
junctions one can employ "orthodox theory" jfii which 
treats the tunnelling in lowest order perturbation theory 
(golden rule). In the strong coupling regime, however, 
this theory collapse since then the transport is not domi- 
nated by sequential tunnelling. While effects of coherent 
tunnelling have been extensively studiedj&^iiLILiSi the 
question of a general description which includes correla- 
tion effects and is valid both for strong and weak coupling 
regimes remains open. 

In this paper, an analytical formula for the current 
through a single level quantum dot (QD) is derived for 
arbitrary on-site correlation strength, beyond standard 
mean-field theory (e.g. beyond self-consistent Hartree- 
Fock, or Hubbard I, approximation (HIA) 13 ' 14 ' 15 ). This 
formula include local correlation effects which pro- 
vide the well known renormalisation of the localised 
leveli 1 ^i 1 ?i 1 Ti 1 ?i 1 ?i 2 9i 21 The theory and formula presented 
are shown to be consistent with previous model results 
for interactions between localised and delocalised elec- 
trons, e.g. atomic limit, non-interacting limit (on-site 
Coulomb repulsion U — » 0), and the strongly correlated 
limit (U — » oo). While these limits have been treated 
several times, and the last limit has been treated in 
detail concerning Kondo physicsj 22 ' 23 ' 24 ' 25 ' 26 ' 27 it is im- 
portant that the here presented material includes these 
limits. It is also important to note that the renormal- 
isation of the localised level often discussed in scaling 
theoryjiSiiSiiiiSiiSiSliSi is included into the present for- 
mulation. 

Already in 1987 Larkin and Matvee\s2£ suggested a 
simple formula for current through a single resonant level, 



Lrfl 



(u - e ) 2 + (T/2Y 



Here, T L ' R , T L + T R = T, defines the coupling between 
the resonant level £o and the leads, whereas /l/_r(w) = 
f(oJ—fJ>L/R) is the Fermi function for the left/right (L/R) 
lead at the chemical potential Hl/r- The formula pro- 
vides a very simple description of the single level QD 
with no on-site interaction, i.e. U — > 0. Since then, this 
description has been further generalised to arbitrary in- 
teractions in the interacting region^SiiS however, given 
in terms of non-equilibrium Green functions of the lo- 
calised levels. Although this formulation provides a gen- 
eral framework for mesoscopic quantum systems in non- 
equilibrium, it nevertheless lacks the transparency found 
in simple single resonant case. It is therefore motivated 
to derive an analytical formula for the current through a 
single level QD for arbitrary on-site correlation U, which 
is not explicitly given in terms of non-equilibrium GFs. 
Thereto, it is desired that this formula goes beyond any 
standard mean-field theory, e.g. HIA, and include effects 
of electron correlations of the localised states. However, 
one should not expect that, for instance, the Kondo effect 
may be included into this formula, since the simplicity of 
the intended expression would be lost by an adequate 
treatment of non-equilibrium Kondo physics. Effects re- 
lated to Kondo physics are hence omitted in this paper. 



The outline of the paper is as follows. The general 
model used for calculations of the current is introduced 
in Sec. |HJ whereas the specific model used in this paper 
is discussed in Sec. IIIII In Sec. IIVI this model is dis- 
cussed for the case of Hubbard's approximation. Then, 
in Sec. [V] the formulation of the problem is discussed 
in terms of many-body operators, where the HIA is re- 
derived as well as a more advanced approximation of the 
system (loop correction), in which effects from electron 
correlations are included. Here, also the formula for the 
current through the system is derived, both in the HIA 
and in the loop correction. Charge and current conser- 
vation of the discussed approximations is proved in Sec. 
fVll a few numerical examples are considered in Sec. IVIII 
and the paper is finally summarised in Sec. IVIIII 
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II. MODEL 



III. SINGLE LEVEL QUANTUM DOT 



In many cases of modelling mesoscopic systems for 
transport, it is reasonable to regard the leads and in- 
teracting region as separate subsystems which interact 
via a tunnelling Hamiltonian Ht- Single or coupled QDs 
attached to leads are examples of often studied systems 
where such an approximation is appropriate, where the 
generic model is given by 



H = TLl + Hr + H 



int 



Ht- 



(1) 



Since the properties of the interacting region, 7i; n t, of- 
ten governs the transport properties of the system as 
a whole, the lead Hamiltonians are normally modelled 
as simple non-interacting electron gases, e.g. H-l/r = 

Yjkaeh/R £ ka c \ a \, wnere 4<r ft) creates (annihilates) an 
electron in the left/right (L/R) lead at the energy Ek a 
and spin a =f, j- The tunnelling interaction, Ht, ac- 
counts for the tunnelling of electrons between the leads 
and the interacting region and in the simplest case only 
single electron tunnelling is taken into account. Hence, 



l 'k 



ntr^kcr"""'^ 



H.c), 



(2) 



kner 



where Vka is the hybridisation between the localised and 
de-localised electrons in the leads and interacting region, 
respectively, and d\ l(J (d na ) creates (annihilates) an elec- 
tron at the nth level in the interacting region. 

The model for the interacting region, H m t, may be of 
a more or less complicated structure. Nevertheless, it is 
normally possible to separate it into a single electron zero 
Hamiltonian term and an interacting term, i.e. H- m t = 
H a + Hi, where H = J2, la - £ nad\ ia d na , and where Hi 
includes electron-electron interactions that take place in 
the interacting region. 

With these basic assumptions, along with the anti- 
commutation relations {k, c\, a ,} — 6kk'Sa-a-', {^<t,4'} = 
8 aa i, and all other equal to zero, it is possible to derive 
an exact expression for the stationary current through 
the system^ 



e 

2h 



tr Im J |[r L -T r ]G<{oj) 



+ [f L {u)T L - f R (<j)T H ][G r (u) - G°(w)]dw. (3) 

Here > T nn R <j = 27T T 1 keL/R V *kna V kn'J{^ ~ £ ka) 18 the 

coupling between the leads and the interacting region, 
whereas G <,r,a (a;) are the lesser, retarded, and advanced 
forms of the GF matrix for the interacting region. 

The model for the current given Eq. © is the ba- 
sic equation for the calculations of transport properties 
in this paper. In the rest of the paper, the GF of the 
interacting region will be investigated for single QD in 
different approximations. 



Consider a single level QD with on-site Coulomb re- 
pulsion U , e.g. 



Hi„t = 2J s dld„ + U 



(4) 



where Eq is the single electron level and n a = d\d a . In- 
troduce the GF Ga(t, t') = (-i){Td a (t)dl(t')). The equa- 
tion of motion for the operator d a is given by 

d \ 

i-F^ ~ £oj d a = U rigda + y] v% a \, (5) 

where a is the opposite spin of a. Further, the equations 
for n^d a and k are given by 



0_ 

Of 



Eq-U \n s d a = y^(— Vkadydycl- 
' k 

+v* ka n < j\ + vl s d a dlcks), (6) 



and 



ig^ ~ Ska )k = Vkada, 



(7) 



respectively. These equations lead to a coupled sys- 
tem including both one- and two-electron GFs, like for 
instance (-i)(Tn„(t)d a {t)dl{t r )) , (-i)(Tk(t)4(*'))i etc - 
In order to proceed analytically (and/or numerically), 
one has to resort to one or another approximation. This 
will be done in the next subsection. However, the aim 
of this paper is to establish the equivalence between the 
described method and the method of using many-body 
(Hubbard) operators. Thus, before embarking into the 
details of the approximation schemes, the concept of us- 
ing Hubbard operators will be introduced. 

As is well known, the Fermi operators d a {d\) can be 
expanded in the eigenstates of the interacting region. In 
the present case this is done by introducing the empty, 
singly, and doubly occupied states, e.g. |0), \u), and 
|2) = | Tl)i an d the outer products (projection operators) 
X pq = \p)(q\H By resolution of unity one has 



da = Y / \p){p\d a \q)( q \ = '£(j>\d a \q)X^ 



{0\d a \a)X" a + (vK\2)X 



ir2 



(8) 



and 4 = {a\d a \0)X°° + (2|d CT |a)A 2 *, giving d\da = 
X aa + X 22 and n^nj, = X 22 . Henceforth, Fermi-like 
transitions (changing the total number of electrons by 
an odd integer) will be denoted by X pg , whereas Bose- 
like transitions (changing the total number of electrons 
by an even integer) are denoted by Z pq , while diagonal 
transitions are denoted by hP = Z pp . The single level 
QD can then be written as 



Hi n t — E p h p , 



(9) 
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where E — 0, E a — £q , and E 2 = 2e + U . 

In terms of the Hubbard operators, the tunnelling 
Hamiltonian takes the form 

H T = $>* CT cU* 0<T + VaX 92 } + H.c), (10) 

ka 

where r\ a — { cx | cZ CT 1 2 ) (77^ = ±1) and (0\d a \a) — 1 ac- 
count for the selection rules between the different tran- 
sitions. In a similar way, the GF for the QD can be 
expanded according to 

Qa{t,t') = 

= H)(T(A° CT + Va X s2 ) (t)(X°° + r) a X 2 °)(t'))u 
= G 0(r (t, t') + r\ a G s2aVl {t, t') 

+V*G 0lT29 (t,t') + G 92 (t,t'), (11) 

(the subscript U signifies the dependence of the gen- 
erating functional, to be used later) where Go a = 
Gocro-o and Gg 2 = Gg 22 g. In this example, the 
generating functional is defined by the action S — 

cxp H Ito~' H'{t)dt], 31 i 32 i 33 where the disturbance po- 
tential is given by 

n'(t) = u (t)h° + ^u rTa ,{t)z™' + u 2 (t)h 2 . (12) 

aa' 

The equation of motion for the operators X 0(7 and X a2 
are given by 



d_ 

dt 



1)2 



* U 



+V* ka (h° + h°)\ + V* k fZ™C ki r), 



(13) 



and 



dt 



i^r-AL )X 



0-2 



fc 



t 



+r] a vl a {h s + h 2 )k + i ls v* ks Z s °c ks ), 



(14) 



where A° = E a — Eq and A 2g . — E 2 — E s are the bare 
transition energies in the interacting region. In the fol- 
lowing treatment, terms like c* ks Z 02 are neglected. These 
terms give rise to propagators similar to those considered 
in the theory of superconductivity^^ which is beyond 
the scope of the present paper. 



{~i){Tn s {t)d a (t)dl{t)) is given by, c.f. Eq. ©, 
Uj t ~ £0 - l/W(t, = S(t - 0(»ff(*)> 

k 

k 

k 

The HIA corresponds to the de-couplings 15 

H)(T(n ff k)(i)4(0) = (n<r(t))F kcr (t,t'), 

H)<T(^ CT 4)(i)4(t')> = 0, 

{-i)(T(d a dUk S ){t)d\(t')) = 0, 
where F ka (t,t') = (-i)(Tk(i)4(*')), which yields 

^ - £ o - t/Wfo - S(t - f)<n*(*)> 

Using Eq. Q), the equation for the transfer GF F ka (t, t') 
can be integrated to 



F ka (t,t') = v ka 



t o -i0 



g k *(t,t")SAt",t')dt", (15) 



where g ka (t,t') = (-i)(Tk(t)cl a (t')) is the GF for the 
electrons in the leads satisfying (id/dt — e k(7 )g kcr (t, t 1 ) — 
6(t — t'). Hence, the (Fourier transformed) QD GF 
Ga(t, t') reduces to the well known result 



iw — £q — (1 - {ng))U 



[ilu - e - V a ][iuj -e -U]- {n s )UV a 



, (16) 



where V a = V a (iu) = J2 k \v ka \ 2 /{iu ~ e kG ). 

The HIA for arbitrary U contains three important limit 
results for, namely, the atomic limit v ka — > 0, the non- 
interacting limit J7 — > and the strongly correlated limit 
U — > 00, e.g. 

hm G a {tLu) = - \- -, (17) 

v ka ^0 VjJ — Eq ILU — £q — U 



IV. HUBBARD'S APPROXIMATION 



lim g a {iw) = 77^: (18) 

a—o iuj - e - V a (iuj) 



One of the simplest approximations of the GF for the 
QD is obtained within Hubbard's approximation of the 
two-electron GF, henceforth referred to as the Hubbard 
I approximation (HIA). For completeness, the result is 
derived in both the Fermi operator and the Hubbard 
operator representations, in order to elucidate the cor- 
respondence between the two approaches. The equa- 
tion of motion for the two-electron GF Gaa{t,t') = 



e.g. the result from the exactly solvable Fano- Anderson 
model )? 5 '?^ andi 5 - 

(7^oo lu! — Eq — (1 — (n 5 -))l / cr(«w) 

respectively. The results in Eqs. lO, l[T7|l. lfT%|) . and l(T^|) 

will now serve as a basis to compare the analogous result 
obtained within the approach with Hubbard operators. 
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V. FORMULATION IN TERMS OF 
MANY-BODY OPERATORS 

The expansion of the operator d„ in terms of Hubbard 
operators, Eq. (JSJ, provides the basic rules of the expan- 
sion of the QD GF Q a given in Eq. In the latter 
expansion the propagators are constructed to depend on 
the generating functional, Eq. (|12fl . through the defini- 
tion 

= Jts) ' 

and similarly for the other GFs. In taking the time 
derivative of the GF with respect to t one must also dif- 
ferentiate S, i.e. (c.f. Ref. 132) 



4tSI"(() = T 
dt v ; 



= TS[x 0,T (t),n + n'(t)}. 

With the disturbance potential given in Eq. (|12|l the 
commutator [X a ,H'], a — Oct or ct2, becomes 

[X°\H'} = AU a0 (t)X ° + U arf (t)X°% 
[X 92 ,H'] = AU 2s {t)X s2 - U ar r(t)X° 2 , 

where AU a0 (t) = U aa {t) - U (t) and AU 2s {t) = U 2 (t) - 
Ua{t). Hence, omitting transitions like Z 02 , the equations 
of motion for the GFs Go a a{t,t') and Ga 2 a(t,t'), where 
a denotes the conjugate of any of the transitions Oct and 
ct2, become 



i^-A° a0 -AU a0 (t) )a hT „ (i.l 1 ) 



- U^(t)G om (t,t') = S(t - t')P Qaa (t) 
+ Y,( v ^(-mnih° + hnk(t)X a (t')) u 

k ^ 

+vU-i){T{Z^c ks ){t)X- a {t')) u ^ (20a) 
d .„ 



dt 



AU 2s {t)\G rT2 - a {t,t') 



+ U aS (t)G a2 - a {t,t') = 8{t - t')P s2 - a {t) 

+Y.(wiA-i)mih s + h 2 Mt)x a (t'))i 



+r ls vU-i)(T{Z s °c k ^{t)X a (t')) u 



(20b) 



Here, PoaaW = (T{X° CT , X a }(t)) v and P g2 - a {t) = 
(T{X ct2 , X a }(t))u (end-factors) are spectral weights of 
the respective GFs , playing an important role in this 



formulation of the theory. Also, let Po a = Po a aO and 
P s2 = P&22&, for a shorter notation. Physical quantities 
are drawn out of the involved GFs in the limit U^(t) — > 0, 
since the sources are introduced in order to generate a di- 
agrammatic expansion of the GFs only. 

The present diagrammatic expansion of the GFs is gen- 
erated through functional derivatives of the GFs with re- 
spect to the source fields U^(t) in Eq. (|12|l . where the 
functional differentiation operators arise from the three 
operator propagators (~i)(TZ^(t")k(t)X h (t')) u . To see 
this, consider the variation of the GF, say, Go CT (M') 
with respect to the source fields U^(t), e.g. 5Go a (t,t'). 
Through an analogous procedure as described in Ref. l32l 
one finds that 



{-i){TZt{t")X 0a {t)X a0 {t'))u = 

S 



5U ( (t" 



G 0a (t,t'), (21) 



and likewise for the other propagators. Hence, the equa- 
tions of motion in Eq. I|20|l can be rewritten in terms 
of the normal two operator GFs and functional deriva- 
tives thereof. Using that the transfer GF F k aao{t,t') — 
(—i)(Tk(t)X a0 (t'))u can be integrated similarly as F krT 



giving 



Fkaao(t, t') — Vko 



to-i/3 



g k<J (t,t")[G Qa (t",t') 
+7 h G s2a0 {t" ,t')]dt" . (22) 



and Eq. I|21|) . the equation of motion for the matrix GF 
can be written as 

i^I -A - U(t)) G(t, t') = S(t - i')P(t) 



+[P(t + ) + R(i+)] f° "" 

J t 



V(t,t")G(t",t')dt",(23) 



where / is the identity, A° = diag{A° , A° , A° x , A° T } 
(diagonal matrix) contains the bare transition energies 
E a - E and AS] ff = E 2 - E s . The source 



1 ct0 



fields are contained in U(£) = diag{Ui(£), \J 2 (t)}, where 
U n (t), n = 1,2, arc 2x2 matrices defined by 



Ui(t) = 



_ ( Au ]0 (t) u n {t) 



U u (t) AU i0 (t) 



, . _ ( AU 2l (t) -U n (t) 
2( ) \-Un(t) AU 21 (t) 



(24a) 
(24b) 



The end-factor P(f) = diag{P 1 (i),P 2 (t)}, which arise 
due to the non-commutativity of the Hubbard operators, 
contains the spectral weights of the components. Here, 
each entry P n (t), n = 1,2, is given by 



P 0T (t) PotioW 



PX(*) = o 



01 To 



(*) Podt) r 



p m - ( ^ \ 
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where the vanishing off-diagonal components are inserted 
for completness. It may be noted though, that the off- 
diagonal components in P„, n = 1,2, are non- vanishing 
whenever spin-flip transitions occur in the system. In the 
spin degenerate case the system reduces to a 2 x 2 matrix 
equation, since the spin f and J, equations are equal. 

The functional differentiation operator matrix R = 
diag{Ri,R,2} has been introduced, arising by the same 
arguments as P(t), where R„, n = 1, 2 are defined by 



Ri(t) - 

R 2 (*) = 
with the components 

R0crcr'0{t) = H (\ 
Ra22a'(t) = 1 



Rono(t) Roao{t) 

-Rj.224. (t) Rl22f(t) 
R]'22l(t) R\22\(t) 



6U (t) 8U a , a {t) 



6U ss .{t) 



5U 2 {t) 



(25a) 



(25b) 



Note the order of the spin indices in the second term of 
Eq. (|25a|l . Finally, the tunnelling interaction matrix V 
is given by 



V(M') 



V'(M') ^V'(M') 
a z V'(t,t') V'(f,f) 



(26) 



where V'(t,t') = diag{V T (t, *'), VJ.(t, *')}, with V a (t,t') = 
Sfcei _r \ v ka\ 9ka-(ti t'), and a z is the z-component of the 
Pauli spin vector which accounts for the selection rules 
defined by rj a = (<r\d a \2). 

This concludes the definitions of the equations for the 
QD GF in terms of many-body operators in its gen- 
eral form. The next step is to find valuable approxi- 
mations of the local properties that can be used in the 
non-equilibrium description of the system. 



A. Hubbard I approximation 

The HIA corresponds to omitting all functional 
derivatives, and in addition, putting all the averages 
(TZt(t))u but the diagonal ones, (T(h° + h a ){t))u and 
(T(h a + h 2 ){t))u, equal to zero. Hence, the Fourier 
transformed equation of motion for G in the HIA be- 
comes (as U^(t) — > 0) 

(iul- A°)G(*u) = ~P + PV(iu)G(iu) (27) 
which yields the solution 

P)a 



G a(iu) = 



iu - A° - P 0a V a 



and 



Gct2aoiiu) = Tjo 



P02V0 



iu - A° - P s2 V a 



G 0a (iu)). 



The end- factors Po CT , P s2 , interpreted as spectral 
weights, have to add up to unity, i.e. P$ a + P S 2 = 1 
(proved in Sec. IVI|) . For later use, one also notes that 
the occupation number (n a ) = (h a + h 2 ) = P a2 , hence 
Pocr = 1 — {rig). By simple algebraic manipulations one 
thus finds that 



G 0a (iu}) 



(iu - A° ff - P s2 V a )P 0a 



(28) 



[iu - A° - V a ] [iu - ASy - P 32 UV a 
Therefore, the sum Goo- + r]o-Ga2aO gives the expression 
G 0a (iu) + r] a G s2 ao(iu) = 

(29) 



(iu - A° 29 )P 0<T 



' [iu - A° - V a ][iu - A° 2S ] - p 97 yv a 

Similarly, one finds that the sum rja-Go^g- + G s2 can be 
written as 



ri<?G Qa2s {iu)) + G s2 (iu) = 

(iu - K )p s2 



[iu - A" - V a ] [iu - AO ] - P s2 UV a ' 



(30) 



Thus, one finally arrives at (recalling that A^. = E a 
E = e £mdA 2s = E 2 -E a =e + U), 

G Qa (iu) + ri a G s2cy0 {iuj) + rj a G a02s {iuj) + G s2 (iu) = 

= - A° )(P 0ct + P s2 ) - UP 0a = 

[iu - A° - Va] [iu - A%] - P s2 UVa 

iu-e -(l- (n a ))U 

[iu - eo - Va][iu — eq — U] — (ria)UVa 



(31) 



which is identically equal to the expression for the QD 
GF given in Eq. JTBJ. Hence, the QD GF expanded 
in terms of the Hubbard operator GFs gives exactly the 
same result in the HIA as for the GF given in terms of 
Fermi operators. This implies that the results in the non- 
interacting and strongly correlated limits are recovered, 
as well as the trivial atomic limit on the which the use of 
the Hubbard operators is based. 

In the limit of strong on-site interaction (U — > 00; dou- 
bly occupied state excluded) the result in Eq. (fHl is 
easily obtained in the many-body formulation, since the 
Hamiltonian 7ij n t can then immediately be reduced to 
Hint = J2 P =o a EphP and the tunnelling term to Ht — 
J2k( v kcrX 0a + H.c). In this case it is only necessary to 
solve for the GF Goa, since then G a2 , G a2 ao, G§ a2s = 
which lead to that Q a — Go a . For Go a one directly finds 
that 



G a(iu) = 



Pi 



Oct 



A^-P 0a Va 



that is, exactly the same expression as the one given for 
Qa inEq. (Eg. 

It may be seen from the above analysis that a correct 
treatment of the GFs necessarily leads to self-consistent 
calculations of the quantities G and P involved in Eq. 
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(|27|l (assuming that the effects on V from the QD are 
negligible). This is clear since the end- factors Po a = 
(T(h° + h a )) = N + N a and P s2 = {T{h s + h 2 )) = N s + 
N 2 (U^(t) — > 0), where N p , p — 0,a, 2, are the occupation 
numbers of the corresponding states |0), \a), and |2), 
respectively. These occupation numbers are calculated 
from 



N 
N 2 



±JmJ[G& M - G> (u)]<b>, 
^Im^ J G< 2 (lu)cLo, 



(32a) 

(32b) 
(32c) 



in non-equilibrium, where the lesser/larger form of the 
GFs are defined in the following section. Hence, the end- 
factors and GFs self-consistently depend on one another. 



B. Transport equation in the Hubbard I 
approximation 

Clearly the GF in Eq. I|27[) is given as a Dyson-like 
equation which is better seen by introducing the bare GF 
g satisfying the equation (iojl — A°)g(iw) = P. Then Eq. 
(|27|l can be re-written as 



G = g + gVG. 



(33) 



From this equation it is easy to find the re- 
tarded/advanced and lesser/larger forms of the GF, that 
is 

Qr/a _ gr/a _|_ r/ayr/aQr/c^ 

q</> = G r V</>G a , 

where V< = i[f L {u)T L + f R (uj)r R } and V> = - 
f L {uj)]T L + [1 - f R (ui)]T R }, with T L / R defined such that 
T L +T R = T = -2ImV(w). The expression for G</> 
is found by direct application of the Langreth rules for 
analytical continuation^ to the Dyson equation in Eq. 
(|33|l . The retarded/advanced and lesser GF are then 
inserted into the formula for the current, Eq. J3J), where 
the trace now is taken of the 4x4 matrices, whereas in 
the case of Fermi operator representation the trace only 
runs over the spin indices. To show that the formula 
for the current coincide in the two representations, it is 
useful to introduce G aux , where the superscript stands 
for either of < or r/a. The trace of r Q G aux , a = L,R, 
is then given as 



The last line equals the trace over the spin indices of the 
QD GF in the Fermi operator representation. Hence, the 
formula for the current given in Eq. is valid irrespec- 
tive of whether the Fermi or Hubbard operator represen- 
tation is chosen, as expected. 

Here it is useful to establish a formula for the cur- 
rent that is valid for arbitrary U. From the Dyson 
equation of the GF G it is easily shown that G r 0a + 
vAGhaO + G r 0a2s ] + G r s2 = Ql, c.f. Eq. Cg. Defin- 
ing V; = K - «r CT /2, where A CT = ReVJ and Ta /R = 
27r E fcG L/i? Ha\ 2 5{uJ - e ka ) such that T^ + T R = Y a = 
— 2ImVJ, gives the last term in Eq. J3J as 

tr [f L (u;)T L - / fl (w)r fl ][G» - G»] = 

= -i£[/ £ (uOr£ - fn(u;)T R }T a \g^)\ 2 , (34) 



where 



(W - A° - Foal!)' 



|( W 



A 



vsKu-^-UPnV; 



Likewise the first term of Eq. 
tr [T L -T r ]G<(lu) = 

= *]T(r^ - r«)(/i( w )rj + / fl ( w )r«)|^( w )|^5) 

since tr T a G r (f L T L + f R T R )G a = E a K(h^ + 
f R T R )\G^{uj)\ 2 . Summing the two terms in the current 
amounts to the formula 



J = JE/ T L a T R [f L {u)-f R {u)] 
(oj - A% - P 0a U) 2 



7 dLu. (36) 



|(o, - A° - V')(u - A%) - UP S2 V;\* 

In the non-interacting limit this formula reduces to the 
well known result 

t T^T R 



J 



A 



However, also in the strongly correlated limit, U — > 00, 
the formula for the current becomes particularly simple, 
e.g. 



J 



E 



pipi? p2 



| W -A° -P 0(7 ^P 



[/i(w) - f R (u)]duj, 



since, mathematically the only difference between the QD 
GF Q a in the two limits is the presence of the end-factor 
Poo- in the strongly correlated case. Physically, the ap- 
pearance of the end-factor ensures that the QD is popu- 
lated by at most one electron, as expected. 



tr r Q G aux = 

T^ct /riaux 1 

1 -pa //iaux 

- 1 t v^ot 

-1 1 [Lr 



Gaux 
J.2 

/naux I 
ziaux 
ra^aux 1 riQ^aux 



Gaux ^ 1 pft//^aux riaux \ 
I9tnj i- 1 1 {^ \ —<^t 2 \o) 



)+rj(-G; 



paux 
Ll 0T2J, 



aux 
0J,2T 



Gaux\ 
12 ) 



Gaux\ 
T2 J 



012T 



G 



T2 



C. Renormalisation of the transition energies — 
loop correction 

As was seen in Sec. IV Al the the use of the Hubbard 
operators in the HIA may seem as an undesired compli- 
cation to the problem, due to the presence of the source 
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fields U^(t) and the non-trivial substitution of the three 
operator propagator in terms of functional derivatives of 
the GF. Nevertheless, the HIA served to establish a rela- 
tion between the two different methods and to see that 
they provide equivalent results in the important limits, 
Eqs. fTTjl . fl%|) . and (fl^|l . The analysis presented in this 
section will provide a glimpse of the power enabled within 
the introduced framework. 

From scaling theory one finds that the localised states 
should be renormalisedi^iSiiLiSiiSiSSiSi due to correla- 
tions between the QD states. Within the present ap- 
proach it is possible to obtain the similar result in a 
rather straight forward manner and include this renor- 
malisation of the transition energies into the GFs. The 
result is obtained by effecting the functional differentia- 
tion once to the GFs, neglecting the fluctuations of the 
end-factors. As in the HIA, all averages (TZ^(t))u but 
the diagonal ones are zero. The same hold for the GFs 
Goo-s-o and G a22a although functional derivatives thereof 
may not be zero, as seen below. 

The structure of the equation of motion, Eq. , sug- 
gests that the GF G is given on the form G = DP, where 
D denotes the locator carrying the local on-site prop- 
erties of the GF, e.g. its position and width. Hence, the 
variation of the GF amounts to vary both the locator and 
the end-factor, however, any fluctuation of the spectral 
weight will be omitted which then gives <5G = (<5D)P. 
The locator satisfies the matrix property DD 1 = I = 
D X D, hence ^DD" 1 ) = (JD^D" 1 + D^D" 1 ) = 
which leads to the identity 51) = -D(JD~ 1 )D. Thus, it 
is necessary to study the locator and its inverse. 

In general, the equation of motion for the locator ap- 
pear very similar to Eq. (|23[l . however, replacing G by 
D, and P by / in the first term on the right hand side of 
Eq. . The resulting equation of motion for the loca- 
tor suggests that the inverted locator D _1 can be written 



T>- 1 (t,t')=d- 1 (t,t')-S(t,t'), 

where the bare locator d satisfies the equation [id/dt — 
A°-U(t)]d(t,t') = S(t-t')I, whereas the self- operator 
is identified by 



S{t,t') = {[P(<+)+R(i H 



to-i/3 



V(t,*i)D(ti,t 2 )} 



in 



xT>~ 1 {t 2 ,t')dt 2 dt l . 



In the present case it is sufficient to replace the inverted 
locator by its corresponding bare quantity, i.e. letting 
D 1 — > d 1 giving <5D = — D(<$d _1 )D. This means 
that the diagrammatic expansion of the GF is terminated 
after the first functional differentiation, giving the so- 
called loop correctionMiQSSi^ The above observations 
then lead to SG = — D(<5d)G. Continuing the functional 
differentiation to higher orders generate higher order dia- 
grams that account for additional many-body correlation 
effectsj^iiSi for instance contributions from the Kondo ef- 
fect. 



Application of the functional derivative, say, Roa-a-'o to 
the inverted bare locator gives i?oo-o-'o(£ + )d _1 (i2, £3) = 
—5(t 2 — i3)i?oo-o-'o(£ + )U(£2)- It is then easy to see that 

iWo(i + )U(£ 2 ) = 

= -iS(t + - t 2 ). (37a) 

Scrcr' - ^crT^cr'T -^tr'T^oi 

-6<T'j.<5<rT <W - Su^ui 
Saidc-i S^^Sui 

S a 'lS a 1 



iW'o(£ + )U(t 2 ) 
= i5{t+ -t 2 ). 



(37b) 



5 S a> - fe'i<5<H 

V o o -<M 3 ' T t° 



o 
o 



in general. 

As previously was pointed out, the GFs Goo-s-o = 
and G s22a = 0, which implies that also Go a2a = and 
Ga2ao = 0. Therefore, the functional derivatives Roaao 
and R a22a applied to the GFs give zero contribution, 
whereas Roaao and Rg 22a acting on the GFs give the 
loop correction. Employing the results in Eqs. I|37|l to 
the GF G ai give 

x^D a&0 (t" 1 t+)G Qal (t+ ,f) 
-D a2a (t",t+)G s2l (t+,t')). (38) 

Hence, in the limit Ut(t) — > the Fourier transformed 
components GWa and G s2 a of Eq. i(2"5|) reduce to 



with A s = 
SA a0 



[lUJ - A a0 - P 0a V a )G <ja{iu) = P Qa a 

+[r)< r Pa a V a + ri s 8A 2s ]G s2 a{iuj), 

[lUJ - A 2S - P S2 V a )G s2 a{iuj) = Pa2a 
+ [VaPs2V a + 7/ ( T(5A (T o]Go C ra(«w), 



(39a) 
(39b) 



SA & and 

- E 

>-7T ^ 



2n 



VkS\ 



f{£ks) - /M 

Ska - U 



k£L,R 

x {-2Im[£>5 ff (u>) + V S D r a2s0 ( w )] }dcj, (40a) 



5A 2( j 



~2tt ^ 



\Vks\ / 

J Ska ~ <JJ 

{-2lm[D r 0rf2a (oj) + Vs D r a2 (u;)}}dcj,(40b) 

where D t t{u)) is the retarded form of the locator D a i(iu>). 
The loop corrections to the transition energies arise due 
to kinematic interactions between particles in the differ- 
ent localised states which is induced by the presence of 
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the de-localised conduction electrons. This is a charac- 
teristic feature of systems with interactions between lo- 
calised and de-localised electron statesi ? 1 '?? The effects of 
the loop correction on the transport properties of meso- 
scopic quantum systems (single and double QD, and for 
spin-dependent systems) have been discussed in Refs. 
. Now, Eq. H!ffl[) can be written as 



G 0a (iu) 



(iu - A 2S - P^V^Pao 
Denom CT (ia;) 



- A 2 



where 



Denom CT (iw) = (iu) - A° - V a )(iu - A 2S - SA^) 



-U(P s2 V a - SA a0 ). 



(41) 



Similarly one derives 



(iU - A a0 - P^Vg)Pa2 

Denom (7 (ia;) 



G 0ct2( t(«w) = r ——G s2 {iuj). 

IU) - A a0 - Poa Va 

The result for the QD GF Q a thus becomes (using that 
Poa + P91 = 1, vl = 1. an d V^Vs- = -1) 



definition, Eq. l|4"Ufl. it follows that 



SA a0 + SA 2g 



~ D r QS + i^ s D r a2m - i ls Dl s2a ~ D r a2 
= {(u - A 2a - P a2 VT) 

+Va[VaP<r2Vr + l la SA 90 ] 
-VaiVaPoaVff + Tj^SA^} 

-(u - Ag-Q - P ^)}/Denom;H 
A ff0 - SA S - A 2a + 5A 2a 



Denom^(aj) 



A? 



A° 

Z<7 



U 



DenomJJ.(cj) DenomJJ.(u;) 

Hence, the sum of the renormalisation energies tends to 
zero as U — > 0, as expected. 

In the third limit, U — ► oo, one should not expect that 
the resulting expression of the GF equals the correspond- 
ing result from the HIA, since the HIA does not contain 
the renormalisation of the transition energies. The form 
of the resulting GF, however, should be similar as in the 
HIA, since the energy of the doubly occupied state tends 
to infinity whereas its population number N 2 — > 0. In- 
deed, dividing the GF in Eq. l|4*2|) by U, the numerator 
of the GF becomes 

(iu> - A° - 6A a0 - 5A 2S )/U - P 0a , 

whereas the denominator equals 



Q a (iu>) = 

= G 0a (iui) + r)o-[G a 2o-o(iu) + G Qa 2 & {kS)] + G s2 (iu) = 
iu - Ao-o - SA 2S - UPoa 



Denom CT (iw) 



(42) 



It is necessary to check that the three basic limit 
result are consistent within this approximation. The 
atomic limit is again trivial since the renormalisation of 
the transition energies explicitly depend on the strength 
of the hybridisation v ka between the localised and de- 
localised electrons in the system, c.f. Eq. (|40|) . Hence, 
<5A CT o, SA& 2 — > as v ka — * 0, which then reduces Eq. 
(g2J) to the result from the HIA, e.g. Eq. fTVft . 

The non-interacting limit, U — » 0, is straightforward 
to obtain. The denominator, Eq. (|41|) . Denom cr (iCi;) — > 
(iu - A° - V a )(iu - A° - 5A aQ - 5A 29 ), as U -> 0. 
Hence, in this limit the expression in Eq. (|42|) reduces to 



G a (iw) 



- A° n - SArrn - SA 2a 



(iu 



A 



V a )(iu 



A 



SA 



ctO 



SA 2 



iu - A° - V a 



which is exactly equal to Eq. (JTSJ, as expected. 

Before moving on to the third limit, it relevant to see 
whether the sum <5A ct0 + SA 2S — > as U — > 0, since the 
renormalisation of the transition energies has to be small 
in systems with weakly correlated particles. From the 



(iu - A° a0 - V*)(iu - A° ff - SA a0 - SA 2S )/U 

-{P s2 V a -5A a0 ). 

From Eq. I|40(l it follows that SA a Q and 8A 2a are finite 
for all U and since A 2s jU — ► 1, as U — * oo, the final 
result becomes 

lim g a (iLu) = — P ° CT . 

U^oo IUJ — A ct0 — PouV a 

Formally, this result equals the one obtained in the HIA, 
e.g. Eq. i|T5)l. with the replacement A^ —* A ff0 , and 



A.o = A° 



k£L,R 

j{e ks ) - /(c) . 



-2ln\Dl s (uj)]duj. 



In this limit it is fairly simple to find an analytical ex- 
pression for the population number (n a ), since N 2 — > 
giving (n a ) = P a2 = N a + N 2 ^ N a . From Eq. J32F 
one then obatins (at T = K) 



(n<x) 



Poa 

T7 



— r o >- 



arctan 



a=L,R 



fig - A a0 1 

Po ff r CT /2 2 



This GF (U — > oo) was derived in Ref. |33 and is consis- 
tent with the result in Ref. Letting A^q - * A^ in 
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the GF and population number, shows that these equa- 
tions are consistent with the equilibrium result by Varma 
and Yafet^ as expected. 

As in the HIA, the equations for the GFs and the end- 
factors have to be self-consistently solved for each bias 
voltage in order to ensure an accurate non-equilibrium 
treatment of the system. In addition, the renormalised 
transition energies have to be found from self-consistent 
calculations, since for instance the transition energy A a0 
depend on the all the other transition energies, through 
the dependence of the locators, c.f. Eq. l]4Up. In princi- 
ple, this amounts to defining the GFs and then calculate 
the renormalised transition energies which should be in- 
serted into a redefined GF, from which the occupation 
numbers are calculated. Self-consistency is, hence, re- 
quired for the GF both with respect to the end-factors 
as well as the renormalised transition energies. In this 
sense, the presented solution with the loop correction 
goes far beyond the HIA, since the values of the end- 
factors will be influenced by the renormalised transition 
energies. This is further analysed in Sec. IVIII 



D. Transport equation with the loop correction 

Similarly as in the HIA, the matrix equation for the 
Hubbard operator GFs can be set in a Dyson like equa- 
tion in the same form as Eq. I|33|l , replacing the bare ma- 
trix A by the renormalised matrix A. This implies that 
the same results for the retarded/advanced and lesser GF 
hold also in this case. Hence, the terms in the formula 
for the current, Eq. 0, are again given by Eqs. 134(1 
and (|33|) while now is given by Eq. i|4"2|) . replacing 

iu! and V a by oj and VJ, respectively. Summation of the 
two terms then gives the current expressed as 

J = JE/ Wt/iH-ZflHU^Hi 2 ^. (43) 

In the non-interacting limit this formula reduces to the 
same expression as in the HIA, as expected from the 
discussion above, and in the strongly correlated limit 
(U — > oo ) the current becomes 

p _ f yL T Rp2 



b - A a0 - P 0a Vl 



that is, the same expression as in the HIA apart from 
that the transition energy here is renormalised. 



VI. CHARGE CONSERVATION 

Before demonstrating the charge conservation it is rel- 
evant to show the relation P 0a + P S2 = 1, previously used 
in order to derive Eqs. I|29|l . (|3(J|) . and (|31[1 . A general 
proof for this relation is achieved by a direct calculation, 



i.e. using Eq. (|32|) one obtains 
P()cr + Ps2 = 

= ~i Im ^ / [G " CT ~ °^ + G ^ 2 _ G%2]dU 
= -- Im E f[G r 0cT + G r 92 ]cLu 

= — tr Im J G r {uj)duj = 1, 

since G> G< = G r G Q . 

Any approximate scheme used for transport calcula- 
tions has to be charge conserving, and here it will be 
shown that both the HIA and the loop correction indeed 
are. In the stationary case it is sufficient to check that 
d/dtNp = 0. In order to show this equality, first con- 
sider the temporal derivative of the occupation number 
A^o, whereas the result for the others are obtained in a 
similar way. Hence, 

J^iVo = -i{[h°,n}} = -i{[h°,n T ]) = 

= -2Im^^(^ CT °k) = 2Re5>to*&„o (t,t). 

ka ka 

Using Eq. (1^21 one obtains 



+9^t')[G a oAt\t) + Va G a s2r70 (t',t)} )dt', 



where 



- T^0(±tTt')e~^ t - t '\ 
g< a (t,t') = ifie^e-^-t'l 

This gives the equation of motion for Nq as 



d 



l\Vko 



ka 



x / [f(^)[GUt',t)+V,G% 2a0 (t',t)} 

J — OO \ 

-[G< (T (t',<)+7 ?ff G< 2CT0 (<',i)])e- £ -( t - t ')^' 

fe<7 

x U{e ka ){i2lTa[Gl a {e ka ) + v«G r s2<j0 {e kr7 ))) 
+ [Go*( £ k*) + T]oG< 2a0 (e ka )} 
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Then, using T^ /R = 2n J2keL/R \vka\ 2 5{u - e k<7 ), one 
finally finds that 



±- Yl r«(/ a (a;)(i2Im[GS (T ( W )+j 7CT G; 2(T0 M]) 

+ [G <(a;)+r ?CT G< ff0 ( W )]V 

By the given expressions for the GFs in the loop cor- 
rection, one has that (omitting the denominators for 
brevity) 

~ i(/ £ (w)rj + / fl (w)r*)(w - A 2ff - <5A ff0 ) 
x(w - A ff0 - <SA 2ff - UP 0a )P 0a 

and 

^ImlG^H + ^Gt 2ff0 (a;)] ~ 
- -2r CT (w - A 25 . - <5A ct0 ) 

x(w - A^o - ^A 2ff - UP 0a )P Qa . 

The denominators of the two expressions are equal and, 
therefore, by summing over the left and right contacts 
one finds that 

-§^0 ~ E([/^ r - + /« r *] r - ^ r AfL^ + f R r R }) = o, 



VII. NUMERICAL RESULTS 

In this section some of the qualitative differences be- 
tween the HIA and the loop correction will be analysed. 
The investigation is restricted to a few cases clearly show- 
ing qualitative deviations between the two approxima- 
tions. A complete analysis of such differences, however, 
is beyond the scope of the present paper. 

In both approximations, the GF along with the end- 
factors are self-consistently calculated for each value in 

the parameter space {eo, U, T^ H , V s d, k R T}, where V s d is 
the bias voltage applied over the system. In the loop cor- 
rection, the self-consistent calculations also include the 
renormalisation of the transition energies. The current 
is then calculated from knowing the QD GF, whereas 
the differential conductance (dJ/dV) is calculated as the 
numerical derivative of the calculated current. This ap- 
proach has been chosen since the QD GF depends on the 
bias voltage in a highly non-trivial way 



A. Non-magnetic system 

First, consider the system defined by the model, e.g. 
Eqs. and Q, to be in the non-magnetic limit, i.e. 

L I R 

e a = £o, T a = To/2, and no external magnetic field. 
Then, as discussed, in the non-interacting limit U — > 
there is no difference between the HIA and the loop cor- 
rection. In the opposite limit, U ~ ■+ oo, the difference 
between the HIA and loop correction is the renormalisa- 
tion of the transition energies, i.e. <5A ct0 = A ct0 — A° 
which tend to lower the energy for the localised state 



showing that the occupation number TVo in the loop cor- 
rection is conserved in the stationary regime. This also 
follows for the HIA by removing the renormalisation en- 
ergies in the expressions above. Similar equalities can be 
shown for the other occupation numbers, e.g. N a , -/V 2 . 

In performing the analogous steps for the temporal 
derivatives for N a , A 2 , and by summing up the occu- 
pation numbers one obtains the expression 

j t (N + J2 N a + N 2 ) ~ tr J (tG< (w) 

+ [f L (u)T L + f R (uj)T R ][G r (cu) - G»])da,, 

which resembles the expression for the current, e.g. Eq. 
(0. Hence, using Eqs. and for the traces, re- 
placing T L - T R by T and f L T L - f R T R by f L T L + f R T R , 
one easily shows that the total charge in the loop correc- 
tion, as well as in the HIA, is conserved in the station- 
ary regime. This then shows the current conservation 
Jl = —Jr, where Jl/r, is the current in the left /right 
lead, and the validity of the formula for the current in 
Eq. @ in the given approximations. 




10/3 20/3 10 

bias voltage (eV sd /r ) 



FIG. 1: (Colour online). Transport characteristics of the QD 
for {eo,kBT}/To = {3,0.014} as U — > oo, calculated within 
the HIA (dotted) and the loop correction (solid), a) The 
differential conductance, and b) the current, as functions of 
the bias voltage. 
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5/6 5/2 25/6 (n-f +ni) (bold) of the QD, as functions of the bias volt- 

bias VOltagG (6V sc j/r ) age calculated within the HIA (dotted) and loop correction 

(solid). Parameters the same as in Fig. [5] 



FIG. 2: (Colour online). Transport characteristics of the QD 
for {e , U, k B T}/r = {5/6, 0.5, 0.014}, calculated within the 
HIA (dotted) and the loop correction (solid), a) The differen- 
tial conductance, and b) the current, as functions of the bias 
voltage. The small ripples in the solid line in a) is due to that 
dJ/dV is the numerical derivative of the calculated current. 



|cr). For A° > (A° < 0) it is expected that the cur- 
rent in the low bias voltage regime \V s d\ — > is smaller 
(larger) in the HIA than in the loop correction. This is 
expected since the transition |0)(er| becomes resonant for 
lower (higher) bias voltages in the loop correction, due 
to the renormalisation. Apart from this, the qualitative 
current-voltage (J — V) characteristics within the two 
approximations are expected to be very similar, which is 
verified in Fig. ^ 

In contrast, tuning the system into the regime 1 < 
A° /r < A0 5 /r (or A0 /r < A0 5 /r < -1) and 
0-5 < U/Tq < 1, for low temperatures, one finds signif- 
icant qualitative deviations in the transport character- 
istics between the HIA and loop correction. The cur- 
rent given in the HIA is peaked for bias voltages such 
that either of the chemical potentials 1-J.l/r lies between 
the two transition energies, see Fig. [21 (dashed). Con- 
sequently, there will be a region of a clear negative dif- 
ferential conductance (NDC) between the two conduc- 
tance peaks corresponding to the transitions A^ and 
Ajg.. Such a behaviour of the transport characteristics, 
which is inconsistent with recent experimental data^^i 
is however not expected to occur in spin-degenerate sin- 
gle level systems. Rather, it is expected that there is a 
plateau in the bias voltage range where the one-particle 
state is resonant while the two-particle state is out of 
resonance. This character is captured within the loop 



correction, see Fig. [21 (solid) . 

The qualitative difference of the two approximations 
in this regime may be understood as follows. To be spe- 
cific, consider the case 1 < A" /r < A^/Tq (whereas 
the case A° /r < AS^/To < — 1 can be understood 
from analogous arguments). The renormalisation of the 
transition energies yields A a0 < A° and A 25 - > A^g.. 
Hence, in the loop correction the state \a) begins to pop- 
ulate at lower bias voltages, see N a in Fig. (solid), than 
in the HIA (dashed), whereas the state |2) remains unoc- 
cupied for higher voltages (N2). (Notice, however, that 
the empty and one-particle states has a lower respectively 
higher population almost throughout the whole range of 
bias voltages. This does nonetheless not alter the fol- 
lowing arguments, since it is the variation of the popu- 
lation numbers that give rise to the changes in the re- 
sulting current/differential conductance.) Therefore, the 
one-particle state almost fully saturate for voltages such 
that the two-particle state is out of resonance. In the 
HIA, on the other hand, the one- and two-particle states 
compete about the available population in the QD, since 
the transition energies lie closer to one another. This, 
in turn, leads to a depopulation of the empty state and 
an overpopulation of the one-particle state. In general, a 
high population number in the one-particle state in com- 
bination with a significant reduction of the population 
(No) in the empty state |0), tends to reduce the tun- 
nelling probability of the one-particle stated Hence, the 
overpopulation in \a) affects the transport properties of 
the system such that the current decreases. Eventually, 
for bias voltages sufficiently large so that the two-particle 
state becomes resonant, the population in \a) decreases 
to its "normal" saturation value which leads to an in- 
creasing current. In summary, correlation effects between 
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particles in the different states of the QD tend to remove 
overpopulation of the one-particle state which, in turn, 
removes the NDC for voltages such that \a) and |2) is 
resonant and out of resonance, respectively. 



B. Ferromagnetic leads 

In spin-dependent systems there is a more significant 
difference between the HIA and the loop correction, than 
in the non-magnetic case. This is seen already by study- 
ing the renormalisation, Eq. H4Q(I . from which it is clear 
that the transition energies A CT o, A 2S ., non- linearly de- 
pend on the properties of the spin a subband. Hence, by 
coupling ferromagnetic leads to the QD where, say, spin 
| is in majority, then causes a stronger renormalisation 
of the transitions Ajrj, A 2 j than what is experienced by 
A jo, A2J. That is, the ferromagnetism in the leads in- 
duce a spin split of the transition energies in the QD. 
This fact has been analysed previously, 33 42 for collincar 
and non-collinear alignment of the magnetisation direc- 
tions in the leads. However, these studies focused on the 
large U limit, whereas the present analysis concerns the 
difference in transport properties between the HIA and 
the loop correction for arbitrary U with collinear ferro- 
magnetic leads. 

The effects considered here are related to the effective 
spin-dependent coupling parameters T^ R , These cou- 
pling parameters can be viewed as to account for spin- 
dependent tunnelling probabilities for electrons through 
the left/right interface and/or different density of spin 
t / I electrons in the leads. However, using the effective 
couplings permits a qualitative analysis of the resulting 
transport properties without specifying the actual spin- 
dependence of the leads and/or the tunnelling interfaces. 

The spin-dependence of the coupling T", a — L, R, is 
parametrised in terms of p a = (Tf - T^)/(T^ + Tf) 6 
[-1, 1], letting r« n = r (l±p a )/2, where T = Tf+T". 
No essential physics is lost by this procedure, as was dis- 
cussed by Martinek et al. 26 In terms of the paramaters 
p a , one can study the transport properties of the system 
for parallel (plPr > 0) and anti-parallel (plPr < 0) mag- 
netic alignment of the leads. Parametrizing in this way 
also allows to study the transport properties for sym- 
metric (\pl\ = \Pr\) and asymmetric (jpL ^ pa) spin- 
dependence of the couplings. 

In the present paper, the QD is coupled to a non- 
magnetic and a ferromagnetic lead, specified by p^ = 
and pr ^ 0. As in the non-magnetic case, the bare 
transition energies A" , A^ 3 are spin degenerate. The 
dressed transition energies, however, are spin split due 
to the spin-dependent couplings Y^ R ^^ as can be un- 
derstood from Eq. (|4U|) . and the discussion above. 

Consider the right lead being half-metallic, e.g. pr = 1 
(or pn = —1). Then from Eq. I|40|) . the dressed tran- 
sition energies are then expected to be spin split such 
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that Ajo < Ajo and A 2 | < A 2 j (or Aj < Ajo and 
A 2 j < A 2 j ). While this behaviour of the dressed transi- 



FIG. 4: (Colour online). Transport characteristics of the QD 
coupled to one non- magnetic lead (pl = 0) and one half- 
metallic lead (pa = 1) resulting from the HIA a) and c), and 
the loop correction b) and d). a), b) Population numbers 
N p , p = 0, a, 2, (faint) and (nj + n.4.) (bold), c), d) Differen- 
tial conductance (dJ/dV). Insets show the J — V character- 
istics of the system. Here, {eo, U, fcsT}/ro = {—5, 5, 0.08}. 



tion energies is general, the renormalisation for the tran- 
sition energy A CT o is stronger for large on-site Coulomb 
repulsion 3 ^^ (U — > 00) and weakens as U — > 0. This 
is also the expected property of the renormalisation from 
scaling theory 17 i 18 i 19 i 20 i 21 i 43 The spin split of the dressed 
transition energies then leads to that a current starts to 
flow in the spin j channel for lower bias voltages than in 
the spin j channel. However, in the example presented 
here, the will not be any current in the spin J, channel 
since there are no available spin J, states in the right lead. 
Hence, the total current solely consists of spin } electrons. 
Therefore, one cannot expect any difference in the spin 
polarisation of the current calculated in the HIA and in 
the loop correction. 

Nevertheless, by tuning the system into the regime 
A° /r « 0, A° ff /r„ « (or A%/T « 0, A^/Tq » 0), 
for low temperatures, one finds significant qualitative de- 
viations in the transport characteristics between the two 
approximate schemes. Most interesting is that the cur- 
rent calculated in the HIA shows a resonant peak and 
associated NDC for negative biases (eV s d = hl—^r < 0) 
not far from equilibrium, see Fig. 01 c), whereas this fea- 
ture is completely absent in the current and differential 
conductance resulting from the loop correction, Fig. 0] 
d). In these plots, the peaks around =pl0 and 0, are 



associated with the transition energies At 



TO 



and Ajj, 



re- 
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spectively, where AJ refers to A? in the HIA and to A a 
in the loop correction. 

As in the non-magnetic case, Sec. IVII Al the difference 
in the two results can be traced down to the population 
numbers N p , 0,er, 2, Fig. 0] a) and b). In both approxi- 
mations, the population number has finite values for 
eV st x/T > —10, and tends to zero for eV st x/T < —10. 
This is plausible since the transition energy At lies 
below the chemical potential of the left lead whenever 
eV s d/To > —10, which permits a leakage of spin J, elec- 
trons from the left lead into the QD. Due to this popula- 
tion, there is a finite probability for a double occupation 
of the QD, that is, that the state \2) acquires a popula- 
tion number N 2 > 0. In equilibrium, then, the popula- 
tion numbers JVj — = N 2 in the HIA, which reflects 
the fact that the transition A^ = /1 (= 0). In the loop 
correction iV T > N± « N 2 > 0, since < A 2l < A 2 -\ . As 
the bias voltage is turned on, such that \xi, ~ M_r < 0, one 
expects the spin \ electrons to exit the QD for rather low 
voltages. Hence, the population number N± is expected 
to approach zero, since there are no new spin j, electrons 
entering the QD from the right. Here, the two approxi- 
mations differ in the sense that iVj_ approaches zeros much 
faster in the loop correction than in the HIA, as the am- 
plitude of the bias voltage increases. In addition, as seen 
in Fig. 0]a) and b), the population number N± is smaller 
in the loop correction throughout the whole negative bias 
voltage range, than in the HIA. In turn, the more rapid 
decay and lower values of Ni in the loop correction leads 
to a faster increase of JVf up to its saturation value close 
to unity, as long as A-po not is resonant. Hence, when the 
state I a) is almost fully occupied the current through the 
QD is suppressed since there is almost no weight of the 
state j t) that can transfer electrons from the right to 
the left lead. However, due to the fast decay of N±, the 
current hardly grows larger than its value in the range 
where N-\ is saturated. Therefore, there is no resonant 
peak in the current calculated in the loop correction. 

Apart from the removal of the resonant peak and asso- 
ciated NDC, found in the HIA, there is another distinct 
difference between the resulting transport characteristics 
in the two approximations. For positive bias voltages 
(jiL — Mfl > 0) around 10 mV, the transition A* be- 
comes resonant. Hence, the differential conductance is 
expected to be peaked around this voltage, Fig. |3|c) and 
d), since At being resonant opens a second channel for 
electrons to flow through, in addition to AIm. Opening 
the transition |0)(t | for transport leads to a reduction 
of the population of spin f electrons in the QD. How- 
ever, the population number iVf ~ in almost the whole 
range of positive voltages, see Fig. 0] a) and b). The 
reason is that the one-particle state is fully occupied by 
spin I, thus there cannot be any accumulation of spin f 
electrons in the QD since /ir < A^ < (J,l, meaning that 
any spin j electron entering the QD from the left will 
immediately exit the QD to the right through the tran- 
sition I j}(2|. The reduction in the population of spin f 
electron in the QD is possibly seen implicitly through a 



slight reduction of the population number N%, see Fig. 
H]a) and b). The small redistribution of the population 
numbers of the QD states gives rise to a small peak in 
the differential conductance. Nonetheless, the conduc- 
tance peak is more apparent in the loop correction than 
in the HIA, which again is attributed to the inclusion of 
correlation effects in the former approximation scheme. 

The double peaks in the differential conductance at 
around —10 mV in Fig. |3|c) and d), are due to the sim- 
ilar effects, that is, the transition |0}(er| becoming res- 
onant which leads to a redistribution of the population 
numbers. In addition, the double peaks in the differential 
conductance around zero bias voltage calculated within 
the loop correction are due to the spin split of the tran- 
sition energies, here A 2 i < A24. 

Finally it is worth to notice that other regimes have 
been considered elsewhere, for instance, the empty or- 
bital regime < A° < A^g. for similar arrangements of 
the spin-dependent couplings The result found from 
these studies are well confirmed within the present ap- 
proach (not shown here), as expected, since the present 
theory goes far beyond both any master equation ap- 
proach or, as shown, the HIA. Nevertheless, one notices 
that the transport characteristics calculated within the 
loop corrections in general give slightly lower height of 
the resonant current peak (shallower NDC in the differ- 
ential conductance) than what is obtained in the HIA. 
The arguments for this character are the same as given 
in the analysis of the above examples. Thus it seems as 
correlation effects tend to reduce, or completely remove, 
features like NDC in the current-voltage characteristics 
of the QD system. 



VIII. SUMMARY 

An analytical formula for the current through a single 
level QD was derived beyond standard mean-field theory 
(HIA) for arbitrary on-site correlation strength. Using a 
diagrammatic technique based on the atomic limit prop- 
erties of the interacting region, enabled an expansion of 
the QD GF in the strongly coupled regime. The local 
properties, e.g. QD GF, are solved in a self-consistent 
fashion with respect to the transition energies and on- 
site population numbers. The derived formula is consis- 
tent with previous results for single resonant levelS&SI in 
the non-interacting limit (U — > 0) and with results from 
strongly correlated systems (U — > 00) 

By means of the derived formula it was shown that 
resonant current peaks and associated NDC found in 
the HIA are removed by effects from electron correla- 
tions that are included into the present description (see 
Sec. IVIIf) . In the non-magnetic case, the NDC in the 
HIA is found to arise due to an exaggerated accumula- 
tion of electron density in the QD for bias voltages such 
that one transition is resonant whereas the other is not. 
The exaggerated population, in turn, leads to a reduced 
conductivity of the available transitions. However, the 
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overestimation of the QD electron density is removed in 
the loop correction which then leads to a plateau in the 
J — V characteristics, in agreement with recent experi- 
mental data42*2i In the spin-dependent case, a similar 
overpopulation of one of the spin states \a) in the HIA 
reduces the transmission through the QD. Again, this en- 
hanced population of the QD is removed in the loop cor- 
rection, whereas the resonant current peak in the J — V 
characteristics vanish. 
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